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Starting from a NJL-type model with N fermion species fermion and difermion condensates and 
their associated phase structures are considered at nonzero chemical potential /i and zero tempera- 
ture in spaces with nontrivial topology of the form S^ ® S^ ® S^ and _R^ ® 5^ . Special attention is 
devoted to the generation of the superconducting phase. In particular, for the cases of antiperiodic 
and periodic boundary conditions we have found that the critical curve of the phase transitions 
between the chiral symmetry breaking and superconducting phases as well as the corresponding 
condensates and particle densities strongly oscillate vs A '--^ ^/L, where L is the length of the cir- 
cumference S^ . Moreover, it is shown that at some finite values of L the superconducting phase 
transition is shifted to smaller values both of /j, and particle density in comparison with the case of 
L = oo. 
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I. INTRODUCTION 



Last years, great theoretical efforts are devoted to the understanding of the QCD phase diagram. Since at rather 

^ small values of baryonic density weak coupling perturbative QCD methods are not applicable, usually effective field 

Y^. theories such as the Nambu - Jona-Lasinio type models (NJL) [1] etc are invoked for such kind of investigations 

O i' [2tI1- Evidently, at low temperatures and baryonic densities one deals with the hadronic phase. But at growing 

1^ , baryonic density, due to a condensation of Cooper pairs of two quarks (diquarks), there appears a phase transition 

T^ , • to the so-called color superconducting (CSC) phase of QCD (see, e.g., the reviews [Tf), where the color symmetry is 

spontaneously broken down. In particular, it turned out that NJL-type models are well-suited for the description of 

OsJ ' the chiral symmetry restoring phase transition and low-energy phenomenology of mesons [8] as well as of properties 

> ; of CSC quark matter 0. 

0> , Since the most attractive feature of NJL models is the dynamical breaking of the chiral symmetry in the hadronic 

phase, the additional influence of different external factors on the chiral properties of these models was also stud- 
ied extensively. For example, they were used to investigate dense baryonic matter in the presence of external 
(chromo)magnetic fields Q. In p articular, it was demonstrated on the basis of NJL models in diverse dimensions 
that both external magnetic [10] and chromomagnetic [11] fields induce the chiral symmetry breaking. Moreover, 
chiral symmetry breaking in four-fermion models was studied in weakly curved spaces [121 . |l3| and in spaces with 
nontrivial topology, when one or more space coordinates were compactified |14l4l6| . In addition, the properties of 
finite size normal quark matter droplets in the language of the MIT-bag model were investigated, e.g., in the review 
[17| . Recently, it was also noted that the position of the chiral critical end point of the QCD phase diagram, which 
could be investigated in heavy ion collision experiments, depends essentially on the finite system sizes |18| . 
\^ • There is also some progress in the understanding of the influence of different external factors on the CSC phase 

5^ I transition. In this context it is worth mentioning that an external chromomagnetic field induces the CSC phase 
transition [19] and that an external magnetic field leads to the appearance of new magnetic phases in the three-flavor 
color superconducting quark matter |20|]. Moreover, the effect of spaces with constant curvature on CSC was studied 
in [21| . Note also that the stability of finite size quark matter droplets in the color- flavor locked phase was investigated 
in the framework of a bag model using the so-called multiple expansion method [22l | . 

In the present work we shall use an alternative approach in order to investigate superconductivity in dense cold 
fermionic matter placed in a finite volume. In particular, we shall study the Cooper pairing phenomenon in the 
framework of a NJL model describing the interaction of N fermion species in compactified spaces with nontrivial 
topology. As in QCD, this model ensures in the usual R^-space the chiral symmetry breaking at rather small values of 
the chemical potential /i, whereas at large values of fi there appears a superconducting phase due to the condensation 
of difermions. 

In this context, let us recall the well-known fact that spontaneous symmetry breaking in low dimensional quantum 
field theories may become impossible due to strong quantum fiuctuations of fields [23]. The same is also true for 
systems that occupy a limited space volume. However, as it is clear from physical considerations, the finite size in 
itself may in some situations not forbid the spontaneous symmetry breaking, if the characteristic length of the region 
of space occupied by the system is much greater than the Compton wavelength of the excitations responsible for 
tunneling and restoration of symmetry. (Indeed, one may recall here well known physical phenomena such as the 
superfiuidity of Helium or superconductivity of metals that are observed in samples of finite volume) . This idea has 
been discussed for some scalar field theories as well as for NJL models in a closed Einstein universe, for instance, in 24]. 
Similarly, if quantum fluctuations of fields are suppressed when the number of quantum fields N tends to infinity, 



spontaneous symmetry breaking might even occur in a finite volume. Indeed, suppose that in a finite volume an 
effective potential of an 0(A'^)-symmetrical model has degenerate global minima. Then, in D-dimensional spacetime 
the transition probability from one minimum to another is proportional to exp{—NL^^^) at zero temperature, where 
L is the linear size of the system [25J. It follows from this expression that if L and N are finite, the transition 
probability is nonzero. This circumstance ensures the vanishing of the order parameter and, as a result, might lead to 
a prohibition for spontaneous symmetry breaking in a finite volume. However, if A^ — )• oo the transition probability 
vanishes and the spontaneous symmetry breaking is allowed. 

In the present paper, the consideration of Cooper pair (difermion) condensation in restricted regions of space is 
performed in a toy NJL model in the mean field approximation, i.e. in the leading order of the large A^-expansion 
technique at T = 0. In particular, in order to be sure that the obtained results are stable against quantum fiuctuations, 
the composite difermion field must be a flavor singlet (like the composite fermion-antifermion field) which technically 
leads to an iV-factor in the part of the effective action arising from fermion loops and thus guarantees the application 
of the 1/N expansion. It is just by this reason that we demand here 0{N) flavor symmetry as opposed to the usual 
SU{N) symmetry. ^ 

The paper is organized as follows. For comparision, we flrst derive in Section II the expression for the thermodynamic 
potential of cold dense fermionic matter, described by a NJL-type Lagrangian, for the case of R^ space. It is shown 
here that for some fixed values of coupling constants two phases are allowed, the chiral symmetry breaking phase 
(at /i < /ic ~ 0.3 GeV) and the superconducting one (at fj, > fXc)- In Section III the phase structure of the model 
is investigated in the S^ (^ S^ ^ S^ space with periodic and antiperiodic boundary conditions for fermion fields. We 
have found a rather rich phase structure in the (/i, A ^ l/L)-plane, where L is the lenght of each circumference S^. 
It turns out that the boundary between the chiral symmetry breaking and superconducting phases as well as the 
corresponding condensate values and particle densities strongly oscillate vs A at A — > 0. Finally, in Section IV similar 
considerations were performed in the case of R^ ® S^ space topology, where we have found smoother oscillations of 
both the critical curve and the condensates vs A. 

II. THE CASE OF R'^ SPACE 
A. The model and its thermodynamic potential 

Our investigation is based on an NJL-type model with massless fermions belonging to a fundamental multiplet of 
the 0{N) fiavor group. Its Lagrangian describes the interaction in the fermion-antifermion as well as scalar difermion 
channels: 



*7'V'fc V^,*7Vf , (1) 





where /x is a fermion number chemical potential. As it is noted above, all fermion fields ipk [k — 1,...,N) form 



a fundamental multiplet of 0{N) group. Moreover, each field tpk is a four-component Dirac spinor; ^p^ = Cipl 



and V'fc — V'feC' are charge-conjugated spinors, and C = 17^7*^ is the charge conjugation matrix (the symbol t 
denotes the transposition operation). Clearly, the Lagrangian L is invariant under transformations from the internal 
0{N) group, which is introduced here in order to make it possible to perform all the calculations in the framework 
of the nonperturbative large-A^ expansion method. Physically more interesting is that the model (1) is invariant 
under transformations from an Abelian electric charge U{1) group: ipk -^ expia4>k [k — 1,...,N). In addition, the 
Lagrangian is invariant under the discrete 7^ chiral transformation: V-'fe ^ l^i^k, ipk — > —'i'kl^ {k — 1,...,A"). The 
linearized version of Lagrangian ([!]) that contains auxiliary scalar bosonic fields cr(x), ^{x), A*{x) has the following 
form 

- r 1 N N A* - A - 

£ = {^k [l^id, + lii' -a\^k' ^a^ - ^A*A - _[^C.z7^^,] - _[^,Z75^C]^ (2) 

(Here and in the following summation over repeated indices k = 1,...,A'^ is implied.) Clearly, the Lagrangians ([1]) 
and ([2]) are equivalent, as can be seen by using the Euler-Lagrange equations of motion for scalar bosonic fields cr(a;), 
A(a;), A*(a;), which take the form 

a(x) = -2^(^fcVfe), A(a;) = -2|(V;^z75^fe), ^*{x) ^ ~2^{{l>kil''i^'i). (3) 



^ Note the important difference to the case of QCD-Uke NJL models with SU(N) color symmetry, where the usual 1/A'^ expansion cannot 
bo applied to colored diquarks due to the lack of a corresponding N factor from quark loops. 



One can easily see from ([3]) that the (neutral) field (t{x) is a real quantity, i.e. ((t(x))'I' = a{x) (the superscript symbol 
t denotes the hermitian conjugation), but the (charged) difermion field A(a;) is a complex scalar, so (A(a;))''' = A*{x). 
Clearly, all the fields ([3]) are singlets with respect to the 0{N) group. ^ If the scalar difermion field A(a;) has a nonzero 
ground state expectation value, i.e. {A{x)) ^ 0, the Abelian C/(l) charge symmetry of the model is spontaneously 
broken down. However, if ((t(x)) 7^ then the discrete chiral symmetry of the model is spontaneously broken. 

Let us now study the phase structure of the four-fermion model (1) by starting with the equivalent semi-bosonized 
Lagrangian ^. In the leading order of the large- A^ approximation, the effective action 5off (cr. A, A*) of the considered 
model is expressed by means of the path integral over fermion fields: 



exp(z5eff(a,A,A*)) = / ll[d^i][di;i]e^p[i / Cd'^x), 



where 



5eff((T,A,A*)=- / d 



^a^(.) + ^A(.)A*(.) 



^cff. 



(4) 



The fermion contribution to the effective action, i.e. the term S^s in Q, is given by: 

exp(i5cff) = / nMi^flMV'i] exp(^ / [^^i^+Vfe + V^^^"V'^ - ^fe^V'^ - ^kK*^k 
•' 1=1 -^ 

where we have used the following notations ^ 

D^ ^iYd^±pL-i° -ct{x), K* =iA*{x)-f^, K^iA{x)j^. 



(5) 



(6) 



In the following, it is very convenient to use the Nambu-Gorkov formalism, in which for each fixed fc = 1, ..., iV a pair 
of fermion fields ipk and ip^ are composed into a bispinor ^^ such that 



, C\ 



Furthermore, by introducing the matrix-valued operator 

z-i ^^' -^ 



= ^lY. 



(7) 



(8) 



one can rewrite the gaussian functional integral in ([5]) in terms of ^fc and Z and then evaluate it as follows (clearly, 
in this case [dtpk][dipk] = [di/'f] Wfc] =M*fc]): 



exp(z5eff) = J llld^i] exp 1^ y ^kZ^kd^xl ^ J f[[d^i] ^^^P [^ J ^UYZ)^kd^x\ - det^/^{YZ) = det^^^^Z), 

where the last equality is valid due to the evident relation detF = 1. Then, using the relation (j4]), one obtains the 
expression for the effective action: 



Ses{a,A,A*) = -Jd''x 



-a\x) + —Aix)A*ix) 



N 
i — Indet(Z) 



(9) 



Starting from (j9]), one can define the thermodynamic potential (TDP) il{a, A, A*) of the model H]) in the leading 
order of the large- A'^ expansion (mean field approximation) : 



5< 



'off 



cr, A. A*— const 



= -Nn{a,A,A*) I 



dS 



(10) 



Here we have supposed that the quantities a. A, A* do not depend on coordinates x. Moreover, without loss of 
generality, one can set the arbitrary phase of A equal to zero so that A is now considered as a non-negative real 



^ Note that the A{a;) field is a flavor 0(N) singlet, since the representations of this group are real. 

^ In order to bring the fermion sector of the Lagrangian (|2]l to the expression, given in the square brackets of l(5j, we use the following 



well-known relations: 9* = -d^, C^C''^ = -{Y)\ C^y^C'^ = (r'Y = 7^ 
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FIG. 1: The case of R^ space topology: The behavior 
of the critical value of the chemical potential ^^ vs H at 
G = 30.06 GeV-^ 



FIG. 2: The case of R^ space topology: The gaps a and 
A vs /i, where /ic ~ 296 MeV. 



quantity, i.e. A=|A|. As a consequence, the detZ in © and the TDP (ITU|) are easily calculated. Indeed, using the 
general formulae 



det 



A , B 
A , B 



Aci[-AB + AAA-^B] 



and detO = expTrlnO, one can find for the TDP (|T0|) the expression: 

r., .. CT^ A2 iTr.^lnl) 



(11) 



where D = A^ + j^D'^j^D and the Tr-operation stands for the trace in spinor (s) and four-dimensional coordinate 
(x) spaces, respectively. Transferring in dill) to the momentum space representation for the operator Z?, we have 



^^-^^^-h^w + l^^^ 



(2^ 



In 



A^ +'y^{^ + ^7" - cr)7^(^ - /17° - cr) 



fl2) 



where in the square brackets just the momentum space representation, D, for the operator D appears. In four- 
dimensional spinor space the 4x4 matrix D has two different eigenvalues e±, each being two- fold degenerate: 



e± 



(E^f 



pl = {E±^if + A^-pl 



(13) 



where E = \J cr'^ -{- p^ ■ Since Tr^ InZ? = 2 lne+e_, one can integrate in (J12p over po and obtain the following expression 
for the TDP of dense cold fermion matter (more details of this technique are presented, e.g., in |26|): 



a^ A2 

f}(o- A) = — + — 

^ ' ' AG AH 



(fp 
(2^ 



Q{K'-f) 



E 



A 



E-, 



(14) 



where the Heaviside step-function Q{x) has been inserted in order to regularize the ultraviolet divergent integral, and 
A is a cutoff parameter that is usually taken smaller than 1 GeV, i.e. A < 1 GeV. Since the TDP p^ is symmetric 
with respect to the transformations a — )• —a, and /j, — )■ —/it, we suppose in the following that a > and /i > (recall 
also A > 0) . It is important to note that the quantities E^ and E'^ defined in (|13l) are the energies of fermions and 
antifermions (quasiparticles) in a medium, correspondingly. Clearly, each energy level is infinitely degenerated with 
respect to the direction of the momentum p. Indeed, there are infinitely many quasiparticles with the same energy 
but with different directions of momenta. 



B. Phase structure 



In order to obtain the phase structure of the initial model it is necessary to investigate the behavior of the global 
minimum point (GMP) of the TDP ([T4| in dependence on the chemical potential fi. The coordinates of this point are 
usually called gaps. (The a- and A-coordinatcs of the GMP are the chiral and difermion condensates, respectively.) In 



the model two types of the GMPs are allowed, i) (a 7^ 0, 0) and ii) (0, A 7^ 0). The GMP of the i)-th type corresponds 
to the phase with broken chiral 7^-invariance only (it is a so-called normal phase), whereas the GMP of the type ii) 
corresponds to the superconducting phase, in which Cooper pairing of fermions leads to the spontaneous breaking of 
the f/(l) symmetry. 

Throughout the paper, we use in our numerical calculations the value G — 30.06 GeV^^. Moreover, for illustrations, 
let us take the cutoff in the momentum integral in ([T4|) to be A = 650 MeV. For this choice of parameters, the TDP 
(ITi)) then predicts at /i = and H = (the last is equivalent to the constraint A = 0) a corresponding value of 
the chiral condensate which is characteristic to some low energy phenomenological QCD-like NJL models 7]. In the 
case under consideration the phase structure of the model depends essentially on the value of the coupling constant 
H. Numerical calculations show that a.t H < G and sufficiently low values oi fj, < fic the GMP of the TDP p^ 
corresponds to the normal phase of the model. However, at /z > /Zc the superconducting phase is realized. Ii H > G, 
then for all values of /i, and even for /x = 0, the ground state of the system corresponds to the superconducting phase. 
The behavior of the critical value of the chemical potential fie vs H is depicted in Fig. 1. (At the same time, one 
may consider Fig. 1 as a phase portrait of the model in terms of /i and H. Then, below (above) the critical line fie 
the normal (superconducting) phase of the model is realized.) In all subsequent numerical calculations the coupling 
constant H is fixed by the relation H = 0.55 G. Hence, the set of model parameters in our investigations is the 
following: 

G = 30.06 GeV~^ H = 0.55 G, A = 650 MeV. (15) 

As a result, we have for the set p3)) fie ~ 0.3 GeV. Moreover, the behavior of gaps a and A vs /i in this case is 
presented in Fig. 2. It is clear from this figure that at fi < fie the GMP of the TDP P^ has the form {a, 0) (as a 
consequence, the system is in the normal phase), where a « 0.3 GeV, whereas for fi > fie the superconducting phase, 
corresponding to the GMP of the form (0, A) of the TDP, is realized in the system. Evidently, in the critical point 
fie there is a first order phase transition. 

III. THE CASE OF S^ S^ (g) S^ SPACE TOPOLOGY 

In the present section we generalize the previously obtained results to the case of a space with finite volume. 
Evidently, this is a reasonable task, since all physical effects take place in restricted space regions. For simplicity, let 
us suppose that our system is immersed into a box with equal linear sizes, < x,y, z < L. It is well known that in 
this case the task is equivalent to the consideration of the model in the space of nontrivial S^ (E) S^ (E) S^ topology 
with quantum fields satisfying some boundary conditions of the form (here we again simplify the problem, demanding 
identical boundary conditions for all coordinates): 

tpk{t,x + L,y,z) = e"°'ipk{t,x,y,z), '4)k{t,x,y + L,z) = e''^"ipk{t,x,y, z), ■i(;k{t,x,y, z + L) = e'""^ tp k{t,x,y,z) (16) 

where < a < 2, L is the length of the circumference S^, and now each of the variables x, y, z mean the path along 
it. Below, we shall use only two values of the parameter a: a = for periodic boundary conditions and a = 1 for the 
antiperiodic one. 

As a consequence, to obtain the thermodynamic potential r2LQ(CT, A) of fermions moving in a space with nontriv- 
ial topology S-^ (E) S^ (E) S^, we should replace the integration over each momentum in 1^^ by a summation over 
corresponding discrete momenta pna following the rule: 



f d^p -j^ ex, ex, 00 ^ 

/ 7^;^/(Px,Pa,P2) -> -^ ^ ^ ^ /(PfcQ,P/a,Pma), _p„a = "^ (2n -f tt), n = 0, ±1, ±2, ... (17) 



A. The case of antiperiodic boundary conditions 

Applying the rule (|17l) with a = 1 in the expression (|14p , one immediately obtains the TDP of the system in the 
case of antiperiodic ( "a") boundary conditions: 

2 a2 o^3 oo oo oo 

f7i,(a. A) ^ — + — -—Y,Y.T. 0(^' - P'a - via - pD [SI^M + £1^ : (18) 

i=0 fe=0 1=0 



where A = tt/L, E'^^i = ^ {S^ki ± fJ-f + A^, £,ki = \/pL + pL +Pia + '^'^' and p,a = A(2z + 1) etc. 

The quantities £^^^1 ^^ CSl) are the energies of elementary one-fermion excitations (quasiparticles) in a dense medium 
which occupies now a finite volume and is constrained by antiperiodic boundary conditions. (The signs -/+ correspond 
to the energies of fermion/antifermion quasiparticles.) Evidently, both fermion and antifermion quasiparticle energy 
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FIG. 3: Phase structure in the case of S^ ® 5^ ® S^ 
space topology with antiperiodic boundary conditions. 
The numbers 1, 2 and 3 denote the symmetric phase with 
o" = 0, A = 0, the chirally broken phase with a 7^ 0, A = 
and the superconducting phase with ct = 0, A 7^ 0, corre- 
spondingly. Here A = tt/L, ^lc ~ 0.296 GeV, Aq ~ 0.375 
GeV, Ai ^ 0.195 GeV. 



FIG. 4: The gap cr vs A = tx/L at h = 0.18 GeV in the 
case of S^ ® S^ ® S^ space topology with antiperiodic 
boundary conditions. Here Ao ~ 0.375 GeV, Ai ~ 0.195 
GeV, A2 ~ 0.149 GeV. 



levels can be labeled by a triple of integers, (z, fc, I), where i > fc > ^ > 0. Clearly, in a finite volume the degeneracy 
of quasiparticle energy levels are partially (or even totally) removed. For example, each energy level with quantum 
numbers (0, 0, 0) is non-degenerated, the level (1, 0, 0) is three-fold degenerated etc. Now, for each energy level (i, fc, I) 
let us put into correspondence the integer Niki = (2i -I- 1)^ -t- (2fc -f 1)^ -I- {21 + 1)^ and, as a result, the real quantity 



./TJjN, 



ikl • 



Using these definitions, we construct an infinite set of real scales of the form 
Ao > Ai > A2 > • • • > A„ > • • • , 



(19) 



where each scale A^ coincides with one of the above obtained expressions Xiki and vice versa, each quantity \iki is 
equal to some element of the sequence (fT9l) . So, Aq = ^A^/S — Aqoo ~ 0.375 GeV, Ai = ^A^/ll = Aioo ~ 0.195 
GeV, A2 = 7X719 = Alio ~ 0.149 GeV, A3 = ^/WfTJ = Am = A200 ~ 0.125 GeV, A4 = y^ATSS etc. 

Recall that our aim is to investigate the phase structure of the system with the TDP presented in (fT8|). This means 
that it is necessary to study the behavior of the global minimum point (GMP) of the TDP vs ^ and A (or L). The 
structure of the TDP p^ suggests the following strategy for studying the corresponding GMP. Let us divide the 
plane (/i. A) into an infinite set of strips, parallel to the /i-axis: 

cjo = {(m, A) : a > Ao}, wi = {(/i. A) : Ao > A > Ai}, • 



w„ = {(/i, A) : A„-^i > A > A„}, 



(20) 



Due to the presence of the theta- function in (I18|) . each sum is indeed a finite one there. Furthermore, for the values 
of (/z, A) from the strip ojq the argument of the theta-function is negative and hence the term with sums vanishes 
there. We see that in this case the TDP is reduced to the quantity il^ao = (o'^/4G'-|- A^/4iJ) whose minimum lies at 
the point (cr = 0, A = 0). As a result, all the points of the strip wq correspond to the symmetric phase of the model 
(see Fig. 3). If the point (/i, A) belongs to the strip wi, then only the energy levels with quantum numbers (0, 0, 0) 
contribute to the sum in ([T5|) . so the TDP p^ reduces to the quantity i^Lai, 



n 



Lai 



a 
4G 



A^ 
Iff 



8A3 



(V^ 



3A2 



/(^ 



3A2 



(21) 



The form of the GMP of this function depends essentially on the values of (/^, A), so in the strip wi, as numerical 
calculations show, there are three different phases (see Fig. 3): the symmetric phase 1 corresponding to the GMP of 
the form (tr = 0. A = 0), the chirally broken phase 2 with (cr 7^ 0, A = 0) and superconducting (SC) phase 3 with 
(a = 0, A 7^0). ' 



It might occur that for energy levels with different quantum numbers {i,k,l) there arises the same integer Afjj.;. 
energy levels with quantum numbers (1, 1, 1) and (2, 0, 0) we have Afiii=Af200 = 27 etc. 



For example, for the 
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FIG. 5: The gap A vs A at /i = 0.4 GeV in the case of 
5^ (8 5^ ® S^ space topology with antiperiodic boundary 
conditions. Here a « 0.178 GeV, b « 0.279 GeV and 
other notations are presented in Figs 3, 4. 



FIG. 6: The gap a ys X at fi — 0.4 GeV in the case of 
S^ (g) S^ (E) S^ space topology with antiperiodic boundary 
conditions. The notations are the same as in Figs 3-5. 



To study the phase structure of the model in the strip uj2, it is necessary to take into account in the sums (jlSp 
the contribution from the energy levels (1,0,0) in addition. In the strip wa the energy levels with quantum numbers 
(1, 1,0) should be switched on in addition to the previous ones, and so on. However, in each of the strips ujk with 
k >2 only the phases 2 and 3 might exist, so for all A < Ai we have found in the plane (/u. A) just these two phases, 
the phase with broken chiral symmetry and the superconducting one. In Fig. 3 they are arranged below the line 
A = Ai and divided by a zigzag, or oscillating, critical line. The amplitude of oscillations of this line is rather large for 
values of A near the value Ai. However, when A decreases, the amplitude of oscillations becomes smaller and smaller, 
and this critical line, or the boundary between phases 2 and 3, tends as a whole to the point (fici 0) at A — >■ 0, where 
fic is the critical chemical potential at L = oo. Moreover, it is clear from Fig. 3 that at some values of A the SC phase 
is allowed at even smaller values of fj, (up to values fi ~ 0.2 GeV) than it occurs at L = cx). 

In Fig. 4 the behavior of the gap cr vs A at /i = 0.18 GeV is depicted (at this value of ^ the gap A equals to zero). 
In Figs 5, 6 the gaps A and cr vs A are represented at ^ = 0.4 GeV, correspondingly. As it is clear from Figs 3-5, 
the oscillations both of the critical line and gaps are characteristic features of the model in the finite volume. Clearly, 
these quantities oscillate strongly vs A. One should also take into account that the gaps a and A from Figs 4-6 are 
really discontinuous functions vs A in the points Aq, Ai, A2 etc. 

Finally, let us discuss the influence of a nontrivial topology on the values of particle density n^(A) = —dQ,La/dii in 
the SC phase. In Fig. 7 its behavior vs A is presented at fixed /i = 0.4 GeV in comparison with the particle density at 
A = (L = 00). It is clear from this figure that at some finite values of L which correspond to values of A from rather 
small vicinities of Afc {k — 2,3, ..) the superconducting phase is realized at smaller particle densities than at L = 00 
(in these cases the ratio n^(A)/n^(0) is less than 1) and at the same value of /i = 0.4 GeV. It was mentioned above 
that the SC phase may occur even at /j, < /j,c if -^ is finite (see Fig. 3). It is interesting to note that in these cases the 
particle density can also reach rather small values. For example, it turns out that the point (/i = 0.21, A = 0.1255) 
GeV lies in the SC phase (see Fig. 3). Numerical calculations show that for this set of parameters the density of 
particles 'n^{X) inside SC matter is much smaller than Uc, where Uc is the density at /j, = ^c and L = 00, since in 
this case we have «^(A) « 0.36 ric- These facts might be understood by taking into account the results of the papers 
[iTI. 122(1 ■ where in particular it was shown that the actual baryonic chemical potential, i.e. the energy per one baryon, 
increases with decreasing size L of a system. Hence, for rather small values of the chemical potential /u, for example 
at /I < /ic, by decreasing L, it is in principle possible, to reach the value of the actual baryonic chemical potential at 
which the SC gap A is opened, although in R^-space the superconducting phase is not realized at the same value of 

/i < ^c- 

Note that the above mentioned as well as the following results should be taken with caution when A > A. The 
reason is that in this case the infrared cutoff A comes closer to the ultraviolet one A, and thus the phase space is 
decreased due to the presence of the 6*- function, e.g., in the expression ((T5)) . 
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FIG. 7: Ratio of particle densities ^' .^^ vs A at /x = 0.4 

GeV in the case of S^ (i) S^ ® S^ space topology with 
antiperiodic boundary conditions. (The quantities A^ are 
presented in the text below ea. (fT9)) .') 



FIG. 8: Phase structure in the case of 5^ ® 5^ ® S^ space 
topology with periodic boundary conditions. Here Ai = 
0.325 GeV, A2 ~ 0.223 GeV and other notations are given 
in Fig. 3. 



B. The case of periodic boundary conditions 

Applying the rule pT)) with a = in the expression p4)) , we immediately obtain the TDP of the system in the case 
of periodic ("p") boundary conditions: 



f^Lp(cr,A) 



a 
4G 






A3 



E E E^i^'-pl- 



Pkp 



Pip) 



cP+ 
'-'Aikl 



_ cP~ 

'^Aikl 



(22) 



■i— — 00 A;— — 00 / — — 00 



i£^kl±^i) +A^,£. 



P%+Plp+Pfp + '^^^ andp.p 



2iX etc. In the periodic case it 



where A = tt/L, E^-f.^ = y {'^iki =c A^j + ly-, au 

is very convenient to separate in each of the sums in ([22]) the contributions from zero modes, so the expression ((22 
can be rearranged in the following way: 
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i=l k=l 



cP+ 
'-'AifcO 



fP- 

'-'AifcO 



8A 



1=1 
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'^Aikl 



fP- 

'-'Aikl 



.(23) 



Now, as in the previous section, we interpret the quantities i£'^7fc/ i^Aiki) ^^ ^^^^ expression as the energies of the 
fermion (antifermion) quasiparticles which can be labeled again by a triple of integers (i,k,l) with i > k > I > 0. 
Clearly, the quasiparticles of the type (0, 0, 0) always contribute to the expression (|23p . whereas the contribution of the 
other energy levels depends on the A-values. So, it is convenient, as in the case with antiperiodic boundary conditions, 
to divide the (^, A) plane into strips 



wi = {(^, A) : A > Ai}, UJ2 = {(/i, A) : Ai > A > A2}, 



{(^, A) : Xn-i > A > A„}, 



(24) 



where for the real quantities Afe we use the same notations as in Section [ill AI which now however take other values, 
i.e. Ai = 7X74 = 0.325 GeV, A2 = v^AVS « 0.223 GeV, A3 = ^/^^^ ~ 0.188 GeV etc. ^ The regions uj,, in ^ 
are constructed in such a way that in the strip cji only the fermion and antifermion levels (0,0,0) contribute to the 
TDP ([53]), where it looks like 



nLpii<T,A) = ^ + ^ - ^ [V(^ + m)2 + A2 + V(^- m)2 + A2' 



(25) 



^ As in the antiperiodic case, in the periodic one the real quantities Ajj,; = \/A^N~^, where Nn.i = A{i^ + fc^ + fi), might be associated 
with corresponding energy level {i,k,l). The real expressions Ai > A2 > A3 > ■ ■ ■ in II24I I are just the quantities Xi^i arranged in a 
decreasing order. 
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FIG. 9: Phase structure in the case of R^ ® S^ space 
topology with antiperiodic boundary conditions. Here 
Afe = A/{2k + 1) {k = 0, 1,2,3) and other notations are 
given in Fig. 3. 



FIG. 10: The gap a vs A at ^ = 0.2 GeV in the case 
of R^ ® S^ space topology with antiperiodic boundary 
conditions. Here Xk = A/(2fc + 1) (fc = 1, 2,3). 



(In this region there are no terms of the infinite sums of (j23p which supply a nonzero contributions to the TDP.) In 
the strip uj2 the energy levels (1,0,0) of fermion and antifermion quasiparticles are switched on in addition, in the 
strip LLI3 the energy levels (1, 1, 0) are switched on in addition, etc. Investigating the global minimum point behavior 
of the TDP (|23p in each of the strips (1^^ . it is possible to obtain the phase structure of the model in the whole (/i, A) 
plane, which is presented in Fig. 8. The behavior of the gaps a, A and particle density n^{X) vs A in the periodic 
case are qualitatively the same as in the case with antiperiodic boundary conditions (see Figs 4-7). 



IV. THE CASE OF R^ » S^ SPACE TOPOLOGY 



In the present section we continue the investigation of the difermion condensation in the spaces with nontrivial 
topology, this time when it is of the form R'^ ^ S^. For simplicity, it is supposed here that the z-axis is compactified 
and fermion fields satisfy some boundary conditions of the form (the x,y coordinates are not restricted); 



'ipkit, x,y,z + L) = e*''"Vfc(^, x, y, z). 



(26) 



As in the previous section, we shall use only two values of the parameter a: a — Q for the periodic boundary 
condition and a ~ 1 for the antiperiodic one. Recall that L is the length of the circumference S^ . Note also that 
the consideration of any physical system in the above mentioned space topology is equivalent to a restriction of the 
system inside an infinite layer with thickness L. In this case, to obtain the thermodynamic potential J7l„(M, A) of 
the initial system, one again simply replaces the integration over p^ in (1141) by an infinite series, using the analogous 
rule: 



dp3 

27r 



/(P3) 



n— — 00 



fiPna), 



Pna = j{2n 



a), 



0,±1,±2,, 



(27) 



A. The case of antiperiodic boundary conditions 



Applying the rule (l?f)) with a = 1 in (ITi)) . one has for the TDP ft^a of dense cold matter in a space of R^ x S^ 
topology with antiperiodic boundary conditions the following expression 



O . AN ^' ^' 2A^ 

i=0 



(Pp 
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ca+ 
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(28) 



where A = tt/L, £^j = J {£i ± ^) + IS? , Ei — \J]P- + p|^ + a"^ , and pia = A(2i + 1). RecaU that the quantities f^j 
(£"^"1) are the energies of fermion (antifermion) quasiparticles, which are labeled by a discrete index i (i = 0, 1, 2...) 



and by a continuous quantity \p\ . Integrating in 

rJLa(CT,A) 



u 
AG 



over two-dimensional momenta p, we obtain 

A2 A „ /A-A^ 
Iff 



67r2 I 2A 



$a(a,A), 



(29) 
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where 



$a(a,A) = {Na + 1) <j 3m (VA2 + ct2 + ^) V (^^' + fT2 + ^) + A2 - 2 (^A^T^ + m) +A 



-|3/2' 



V^LT^+Ai) +A2 



-|3/2 



3/iA2 In 
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+ ^2 
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^ + m) 
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+ ^2 


Y' ( VpL + <^ 


+ A2 



3a*(v^LT^ + m) y(VpLT^ + m) +A2 

+ (m ^ -m), (30) 



and Na = [^^^r] (recall, [x] means the integer part of a real number x). 

It is clear from the structure of the TDP (PSI - IPO]) that to obtain the behavior of its global minimum point vs /i and 
A it is very convenient to divide again the plane (/i, A) into an infinite set of strips (j20p . this time with Afe — A/(2fc + 1). 
® So, due to the presence of the 8-function in 1^^, in the region ujq, i.e. at A > Ao = A, the TDP (P5)l - (l^51) has a 
trivial form, i.e. 



^Laicr^A) 



(/i,A)Gwo 4G 



A^ 
4iJ' 



(31) 



whose global minimum point (a = 0, A = 0) corresponds to the symmetric phase. In the region wi , i.e. at Aq > A > Ai , 
the integer Na from p0[) is equal to zero, hence only quasiparticles with the energies £^q contribute to the TDP in 
this case. Studying the behavior of the GMP of the TDP vs fi and A, we conclude that in the strip uji three phases, 
the symmetric phase 1, the phase with broken chiral symmetry 2, and the superconducting phase 3, occur (see the 
part of Fig. 9 that corresponds to Aq > A > Ai). In the strip uj2 we have in the expression (pO| Na = 1, so here the 
quasiparticle energies £^ contribute to the TDP in addition etc. However, our numerical calculations show that in 
each of the strips a;2, ws,--- only two phases, the phase 2 and the SC phase 3, of the model are realized. The boundary 
between the phase with broken chiral symmetry 2 and the superconducting one 3 again oscillates when A — ?> (see 
Fig. 9). It turns out that in addition to the critical curve the gaps a and A oscillate vs A at fixed values of fj, (see Figs 
10, 11). Comparing these figures with Figs 4, 5, we see that the more dimensions are compactified, the more sharp 
oscillations of physical quantities occur. In particular, it is clear from Figs 10, 11 that in the case under consideration 
the gaps a and A are some continuous functions vs A, whereas in the case with S^ ^ S^ (S) S^ space topology these 
quantities have discontinuities in the points Aq, Ai, etc. 

Finally, in Fig. 12 the ratio of particle densities ^'^Lj vs A inside the SC phase is presented at /i = 0.4 GeV. It is 
clear that in small vicinities around Afe (fc = 2, 3, ..) the density at A ^ is less than the particle density at A = 0. It 
seems intuitively clear that the smaller a particle density is, the easier a corresponding state of the system might be 
created. So to reduce the efforts in obtaining the SC phase, one could simply fix the A-parameter not far from one 
of the Afc-values. Hence finite size effects might promote the transition of a physical system into its superconducting 
phase. 

B. The case of periodic boundary conditions 

Clearly, to obtain the TDP of the system il^p in this case, it is necessary to use in (ITil) the rule (P7|) with a = 0. 
As a result, we have 



f^Lp(o-,A) 
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(32) 



where A = tt/L, f^^ = ^ {£n ± fJ^f + A^ , £„ = Jp^+pl^ + a^, p„p = 2nA (n = 0,1,2...). Recah that the 
quantities £^„ i^An) ^^^ ^he energies of fermion (antifermion) quasiparticles, which are labeled by a discrete index n 



One should not become confused by the fact that we use the same notation for the boundaries A^ of these strips in different cases. 
Indeed, as it is clear from the text, the value of each Aj; depends strongly on both the space topology and boundary conditions. 
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FIG. 11: The gap A vs A at /i = 0.4 GeV in the case 
of R^ S^ space topology with antiperiodic boundary 
conditions. Here Afc = A/(2k + 1) (fc = 1, 2, 3). 



FIG. 12: Ratio of particle densities "^ ,,,, vs A at /i = 0.4 

GeV in the case of R^ ® S^ space topology with an- 
tiperiodic boundary conditions. Here Afe = A/(2fc + 1) 
(fc=l,2,3). 



{n = 0, 1, 2...) and by a continuous quantity \p\, in addition. We find it convenient to separate in (|32|) the contribution 
from the zero modes, i.e. the contribution from quasiparticles with n = 0. Integrating in ([32]) over two-dimensional 
momenta p, we have 



f^Lp(cr,A) 



where 
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(33) 



(34) 



(35) 



and Np = [A-1 is the integer part of the real number in the square bracket. 

As in section UlIBl we again divide the (/i, A)-plane into an infinite set of strips ojk p4l) . where in the present case 
we have A^ — A/(2fc) (fc = 1, 2, ...). Then, in the region wi the zero energy levels f^g of quasiparticles contribute to 
the TDP dniD, where it looks like 



f^Lp(CT,A) 



A 



(M,A)e<^i 4G 4iJ 127r: 



:F(a,A) 



(36) 



(here F((T, A) is given in (1341) '). In the region a;2 the next energy levels £^^ are switched on in addition, so in the 
expression for the function $p(cr, A) ([55]) we have Np — 1, etc. It turns out that with each region w^ [k > 2) the 
value Np = fc — 1 in (I35p is associated. As a result, in cuk all the quasiparticle energy levels with quantum numbers 
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FIG. 13: Phase structure in the case of R^ (8 5*^ space topology with periodic boundary conditions. Here A^ = A/(2fc) 
(fc = 1, 2, 3), and the numbers 2,3 denote the chirally broken and SC phases. 

n = 0,l,..,,fc— 1 contribute to the expression for the TDP. Studying numerically step-by-step the behavior of the 
GMP of the TDP (p3)) in the strips ujq, uji, UJ2,---, we obtain the (/i, A) phase portrait of the initial NJL model which 
is presented in Fig. 13. In contrast to the phase portrait in the antiperiodic case (see Fig. 9), it has only two phases, 
the chirally broken phase 2 and the superconducting one 3. 

Notice that in the case with periodic boundary conditions the behavior of the gaps a, A and particle density n^(X) 
are qualitatively the same as in the antiperiodic case (see Figs 10 - 12). 

V. SUMMARY AND DISCUSSION 

In the present paper we have investigated the influence of finite-size effects on the superconductivity (SC) phe- 
nomenon which might exist in dense cold fermionic matter. Note that in the SC phase the U{1) charge group is 
spontaneously broken down due to Cooper pair (difermion) condensation. It is well-known that quantum fluctuations 
of fields can destroy symmetry breaking in a finite volume, so studying this effect requires to be sure that such fluc- 
tuations and their corresponding next to leading order corrections will not spoil the leading order mean field results. 
Concerning the application of usual QCD-like NJL models to the description of CSC in the mean field approximation, 
a small perturbative expansion parameter guaranteeing the suppression of quantum fiuctuations is absent. Thus, in 
this case there is no confidence that quantum fiuctuations are in general negligible and, particularly in the case of finite 
systems, that they could not destroy color symmetry breaking. These obstacles were the main reason why we decided 
to consider instead of a QCD-like quark model the "toy" NJL model (1) for fermions which seems to us technically 
more adequate for studying the influence of finite volume effects on the Cooper pairing. In fact, in this model at large 
N a small expansion parameter 1/iV appears, so the next to leading order corrections in 1/N are certainly negligible. 
By this reason, quantum fiuctuations of fields become suppressed and cannot destroy spontaneous symmetry breaking 
for the considered fermion system in a finite volume (see also the discussion in Introduction). 

Let us summarize in more detail some of the main results. First, it was shown in the leading order of the 1/A^- 
expansion that at T = 0, L = oo (L is the linear size of the system) there is a phase transition in the considered 
fermion model from the chiral symmetry breaking phase 2 to a superconducting one 3 at the critical value of the 
chemical potential /x = /Xc « 0.3 GeV. Next, we have studied in the leading order over 1/A^ the phase structure of 
this model both in the space with topology S^ ® S^ ^ S^ and R^ (X) S*^ taking into account periodic and untiperiodic 
boundary conditions for fermion fields. It turns out that for all cases the critical line between phase 2 and 3 as well 
as the gaps a, A and the particle density ^^^(A) are oscillating functions vs A ~ 1/L. Generally we found that the 
more spatial dimensions are compactified, the stronger oscillations occur. 

Secondly, it is interesting to note that at finite L the superconducting phase 3 might be realized at sufficiently 
smaller values of the chemical potential /i and particle density n^(A), than ai L = oo. Indeed, as it is clear from, 
e.g.. Figs. 3, 13 for some values of A the phase 3 occurs at ^ = 0.2 GeV or even smaller values. Moreover, as it is 
clear from the discussion at the end of section IIII AI for such sufficiently small values of /i the particle density inside 
the SC phase might be as small as « 0.4 Uc, where ric is the density at which SC is realized ai L ~ oo. Hence, the 
compactification procedure itself might shift the SC phase transition to smaller particle densities. 

Hopefully, the above investigations could motivate other studies of finite size effects and eventually find some 
physical applications. 
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